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ABSTRACT 

We report on a multi-site photonietric campaign on the high-amplitude 5 Scuti star 
V2367Cyg in order to determine the pulsation modes. We also used high-dispcrsion 
spectroscopy to estimate the stellar parameters and projected rotational velocity. Time 
series multicolour photometry was obtained during a 98-d interval from five different 
sites. These data were used together with model atmospheres and non-adiabatic pul- 
sation models to identify the spherical harmonic degree of the three independent fre- 
quencies of highest amplitude as well as the first two harmonics of the dominant mode. 
This was accomplished by matching the observed relative light amplitudes and phases 
in different wavebands with those computed by the models. In general, our results 
support the assumed mode identifications in a previous analysis of Kepler data. 

Key words: stars: individual: V2367 Cyg - stars: oscillations - stars: variables: 5 Scuti 



1 INTRODUCTION 

The high-amplitude delta Scuti (HADS) stars are defined 
as 5 Set stars with peak-to-peak light amplitudes in excess 
of about 0.3 mag. This amplitude limit is an arbitrary one 
and there is no other physically identifiably distinction be- 
tween HADS and other 5 Set stars. The HADS stars rotate 
more slowly than other 5 Set stars with projected rotation- 
ally velocities vsini < 30kms^^. They pulsate in one or 
two dominant modes which are generally assumed to be ra- 
dial modes, but when observed with greater precision from 
space many more modes are visible. A large number of com- 



bination frequencies and harmonies involving the modes of 
highest amplitude are usually present. The pulsational be- 
haviour of HADS stars seems to be interme diate betwee n 
that of multimode 5 Scuti stars and Cepheids (|Bregeill2007l b 
In fact, first-overtone classical Cepheids and HADS stars fol- 
low th e same period - lumin osity relation with no disconti- 
nuity (|Soszvnski et al.ll2008l '). 

The reason why 5 Set stars pulsate in many nonra- 
dial modes while the more massive Cepheids pulsate only 
in radial modes is due to the behaviour of nonradial pul- 
sation modes in the stellar interior. In the deep interior 
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of an evolved giant star such as a HADS star or low-mass 
Cepheid, nonradial p modes behave like high-order g modes. 
As the stellar mass increases, these modes develop an in- 
creasingly large number of spatial oscillations which leads 
to increased radiative damping. Therefore we find fewer and 
fewer nonradial modes in HADS stars of increasing mass 
until only r adial modes are pr esent in what we then call 
a Cepheid (|Dziembowskil [T977I ) . Indeed, observations sug- 
gest that there is a group of HADS which pulsate just in 
the fundam ental and first over tone modes with no nonra- 
dial modes jPoretti et ahllioosl ) (they might be called low- 
mass Cepheids). While this seems to be the main reason why 
Cepheids do not pulsate in nonradial modes, it may not be 
the only reason and the pro blem needs further investigation 
(jMulet-Marauis et al]l2007l ). 

The study of HADS stars is important because it is 
an interesting laboratory for investigating mode interaction 
and nonlinear behaviour as exemplified by the presence of 
combination frequencies. The use of combination frequen- 
cies for mode identification is well- known in pure g-mode 
pulsator s such as the ZZ Ceti sta r s (IWul I2OOII : lYeates et al.l 
[2005. : .Montgomervll2005l '). ISalonal l|2012l) has recently inves- 
tigated the problem in pulsating main sequence stars. While 
in ZZ Ceti stars the combination frequencies arise from in- 
teractions in a convective zone in which certain approxima- 
tions are valid, this is not the case in HADS and it turns 
out that combination frequencies cannot be used for mode 
identification. 

Mode identification is crucial for extracting stellar pa- 
rameters from the pulsations. The assumption that the 
modes of highest amplitude in HADS are radial is based on 
the fact that when two dominant modes are present, their 
period ratio often agrees with the expected period ratio for 
first overtone to fundamental radial mode, P i/Pp, which is i n 
the range 0.77-0.78. The Petersen diagram (|Petersenlll973 ). 
which is a plot of Pi/Po as a function of Pq, is very use- 
ful for diagnostic purposes. This ratio varies du e to rotation 
with chemical composition (jPoretti et al. 2005|) and by cou- 
pling with the nearby quadrupole mode ( Suarez et al.ll200"^ . 
I2OO7I ). The calculated pe riod ratio is also sensitive to the 
opacity (|Lenz et allbOOSl ). 

V2 367Cyg (KIC 94086 94) was discovered in a ROTSE 
su rvey (lAkerlof et al.ll2000D and confirmed as a HADS star 
bv lJin et all (|2003l ') and lPigulski erakl (|2009l '). The star has 
be en observed by the K epler satellite and recently analyzed 
by iBalona et all (|2012l '). The light curve is dominated by a 
single frequency /i = 5.6611 d~^ of high amplitude. Two 
other independent modes with /2 — 7.1490 d~^ and /s = 
7.7756 d~^ have amplitudes an order of magnitude smaller 
than /i. Nearly all the light variation is due to these three 
modes, harmonics of /i and combination frequencies with 
/i, although several hundred other frequencies of very low 
amplitude are also present in the Kepler photometry. 

The star has a much higher projected rotationa l veloc- 
ity {vsini = 100 ± 10kms"\ iBalona et all (|2012l )') than 
other HADS stars, which makes it unusual. Such a large ro- 
tational velocity complicates the analysis and increases the 
likelihood of mode interaction. The ratio /1//2 ~ 0.792 is 
within the expected ratio of second- to first-overtone radia l 
periods which is in the range 0.79-0.80. IBalona et all (|2012l ) 
attempted to model the observed frequencies as radial modes 
without mode interaction, but were not successful. A model 



Table 1. Information on the photometric multisite campaign. 
EUO: Ege University Observatory; BNAO-Rozhen: Bulgarian Na- 
tional Astronomical Observatory-Rozhon; OAVdA: The Astro- 
nomical Observatory of the Autonomous Region of the Aosta Val- 
ley; SPM: Observatorio de San Pedro Martir; RBO; Red Butte 
Observatory. Telescope apertures are in m. The total number of 
nights observed and the start end ending Julian dates of the ob- 
servations relative to JD 24455000 are given. 



Observatory 


Site 


Tel. 


CCD 


BNAO-Rozhen 


Bulgaria 


0.60 


FLI PL9000 


OAVdA 


Italy 


0.81 


FLI PL3041-1-BB 


SPM 


Mexico 


0.84 


FLI Fairchild F3041 


EUO 


Turkey 


0.40 


Apogee Alta-U47 


RBO 


USA 


0.60 


Apogee Alta E47-UV 


Observatory 


Filters 


Dates 


Nights 


BNAO-Rozhen 


UBVRI 


2 


736.34-737.54 


OAVdA 


BVRI 


3 


737.37-739.56 


SPM 


UBVRcIc 


6 


718.76-732.98 


EUO 


BVRI 


17 


736.29-816.52 


RBO 


UBVRI 


1 


741.68-741.93 



with mode interaction was found in which /i is identified as 
the fundamental radial mode and /2 the first overtone radial 
mode coupled with a nearby quadrupole mode. 

Considering the complex nature of the pulsation, the 
question arises as to whether the modes /i and /2 are indeed 
radial and coupled modes as sugg e sted b y the models. The 
model proposed by IBalona et al.l (|2012h would be greatly 
strengthened if independent mode identification could be 
made. Since the star is sufficiently bright and of high am- 
plitude, multicolour photometry can be used to provide this 
verification. The purpose of this paper is to describe a multi- 
site photometric and spectroscopic campaign on V2367 Cyg 
and subsequent analysis of mode identification. 



2 SPECTROSCOPY 

The spectra were used to estimate the stellar atmospheric 
parameters and the projected rotational velocity, v sin i. 
The spectroscopic observations were carried out with two 
instruments. The Coude spectra were obtained with the 2-m 
RCC telescope of the Bulgarian National Astronomical Ob- 
servatory, Rozhen. We observed the star during the night of 
2011 July 27 in three spectral regions: 4770-4970 A, 4400- 
4600 A, and 5090-5290 A. These three regions contain the 
H^ line, MgH A 4481 A and a region rich in Fe lines. The 
instrument uses a Photometries AT200 camera with a SITe 
SI003AB 1024 x 1024 CCD, (24 ^im pixels) giving a typical 
resolution R — 16 000 and S/N ratio of about 40. The expo- 
sure time was 1200 s. The instrumental profile was checked 
by using the arc spectrum giving a full width at half max- 
imum (FWHM) of about 0.4 A. IRAF standard procedures 
were used for bias subtracting, flat-fielding and wavelength 
calibration. The flnal spectra were corrected to the heliocen- 
tric wavelengths. 

We also obtained spectra with the WIRO longslit spec- 
trograph which uses an e2V 2048^ CCD as the detector. An 
1800 lines mm~^ grating in first order yielded a spectral reso- 
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with the results given bv lBalona et al.l (|2012l ) who find Tcff = 
7300 ± 150 K, loRfl = 3.5 ± 0.1, usini = 100 ± lOkms^V 
lUvtterhoeven et al.l (|201ll ) list the effective temperature as 
TcEF = 6810±130 K, which is in disagreement with the above 
values. 



4 PHOTOMETRY 
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Figure 1. The H/3 line fitted with a model with T^g = 7250 K, 
log 3 = 3.5 (cgs) (short dashes). Two other models with T^ff = 
7000 K (long dashes) and T^ff = 7500 K (dots) are given for 
eomparison. 



lution of 1.5 A near 5800 A with a 1.2" x 100"slit. The spec- 
tral coverage was 5250-6750 A. Individual exposure times 
were 600 s. Reductions followed standard longslit techniques. 
Multiple exposures were combined yielding final S/N ratios 
typically in excess of 60 near 5800 A. Final spectra were 
Doppler corrected to the heliocentric frame. Each spectrum 
was then shifted by a small additional amount in velocity 
so that the Na I D AA 5890, 5996 lines were registered with 
the mean Nal line wavelength across the ensemble of ob- 
servations. This zero-point correction to each observation is 
needed to account for effects of image wander in the disper- 
sion direction when the stellar FWHM of the point spread 
function was appreciably less than the slit width. Because of 
these inevitable slit-placement effects on the resulting wave- 
length solutions at the level of <10kms~^, radial velocity 
standards were not routinely taken. 



3 ATMOSPHERIC PARAMETERS 

Model atmospheres were calculat ed using the ATLAS 12 code. 
The VALD ato mic line databa se (|Kupka et al.|[l999i ). which 
uses data from iKurucjj 19931 ). was used to create a line list 
for the synthetic spectrum. S ynthetic spectra were gener- 
ated using the SYNSPEC code (|Hubenv et all 1 1994 fKrtickal 
1 19981 ). We adopted a microturbulence of 2kms~^. The com- 
puted spectrum was convolved with the instrumental pro- 
file (Gaussian of 0.4 A FWHM for the Coude spectra and 
1.5 A FWHM for the WIRO spectra respectively) and rota- 
tionally broadened to fit the observed spectrum under the 
assumption of uniform rotation. 

The best fit to the H/3 and Hq lines was obtained for 
effective temperature Tctt = 7250 K, and surface gravity 
log 5 = 3.5. We used the MgH A 4481 A hue to determine 
the projected rotational velocity. The best match was ob- 
tained with V sin i = 90 kms~^. In Fig. [T]we show the best 
fit to H/9 as well as the fit using two other models with differ- 
ent temperatures for comparison. Our results for the atmo- 
spheric parameters and v sin i are in very good agreement 



The photometric data were used to determine the 
amplitude-colour relation and hence the mode for each fre- 
quency.The multi-site photometric campaign on V2367 Cyg 
was carried out during 2011 June to 2011 September (Julian 
dates 2455726.69 - 2455816.52) in Europe and North Amer- 
ica with five different sites. Although time allocation for 10 d 
at KISO Observatory, Japan was allocated, no photometric 
data could be obtained due to bad weather. The time of 
these observations coincide with ifepfer quarter Q9.3-Q10.3. 
There are over 127780 short-cadence (1 min) Kepler data 
points in this time range. Table[T] summarizes the sites, tele- 
scopes and additional information regarding the campaign. 

Since the aim of this campaign is mode identification, 
multi-colour UBVRI filters were chosen. However, EUO and 
OAVdA do not have the U filter. The U filter, when it is 
used, leads to high signal-to-noise owing to the poor sensi- 
tivity of the CCDs in this wavelength and is therefore not 
reliable. In any case only four nights of U photometry are 
available which is not sufficient to resolve the frequencies. 
Although we show the U observations in mode identifica- 
tion, they should be given very low weight or ignored. All 
R and I filters are of the Johnson/Bessel type except those 
at SPM (Mexico) which are of the Cousins/Kron type. The 
difference is of little consequence as it has a negligible effect 
on the amplitude ratio and phase difference which are used 
for mode identification. 

When the campaign was initiated in 2011 June, all sites 
observed the target for 3-4 nights. Subsequently, only EUO 
continued to observe the target in July, August and early 
September. A total of 181 hours were spent on observing 
V2367 Cyg within the 98-d duration of the campaign. Data 
were reduced using standard IRAF routines including correc- 
tion for dark, bias an d flat-fielding o f each CCD frame. The 
DAOPHOT II package jStetsonl [19871 ) was used for aperture 
photometry. The observing times were converted to helio- 
centric Julian Date (HJD). 

Since Kepler photometry was available, it was easy to 
check the variability of nearby stars. We used the compar- 
ison stars GSC3556 1047 (CI), CSC 3556 1343 (C2) and 
CSC 3556 627 (C3) because they showed least variability in 
the Kepler photometry (not more than 50 ppm). Because 
not all instruments have the same fleld of view, C2 was not 
observed by some sites and CI was chosen as the primary 
comparison star. The data was therefore reduced relative to 
CI. Outliers produced by poor seeing or transparency vari- 
ations were removed by visual inspection. 

The passbands of corresponding filters in the instru- 
ments on different sites are not exactly the same and neither 
are the optical characteristics of the instruments. Therefore 
there are invariably zero-point differences in each filter for 
each site. Fortunately, we know the frequency composition of 
the vari ations in V2367 Cyg from the analysis of the Kepler 
data bv iBalonZ et all (|2012i ). We can assume that the light 
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Figure 2. Bottom panel: periodogram of the ground-based V 
observations of V2367 Cyg. Top panel: periodogram of the Kepler 
observations over the same time interval. 



variations are well modeled by a Fourier series comprising 
these frequencies where the zero point will vary from site to 
site. This is easy to do by fitting the following function to 
the data by least squares: 



Vi + ei^ 51/ ^'J'^-'' 51/ sin(27r/fcfi 
j=i fc=i 



- <t>k) 



(1) 



where t/i is the magnitude at time ti and is the unknown 
observational error. The number of different sites is s and 
5ij = unless the i-th observation is from site j, in which 



case 5i 



1. The last term is a truncated Fourier series with 



m frequencies which represents the variability of V2367 Cyg. 
By minimizing ^ ef for 1 ^ i ^ N {N being the total num- 
ber of observations), the individual zero points for each site, 
aj , and the amplitudes and phases of the Fourier series can 
be found. 

In Fig. [2] we show the periodogram of the V observa- 
tions and also of the Kepler observations over the same time 
interval (JD 2455718.0-2455817.0). The huge advantage of 
an uninterrupted time series is very clear in this figure. 

We used the first three independent frequencies of 
highest ampl itude obtaine d by the frequency extraction 
described in iBalona et al] (|2012l ): /i = 5.661058, /2 = 
7.148949 and /a = 7.775566 d"^ together with the first three 
harmonics of /i (2/i , 3/i , 4/i ) for the least squares Fourier 
fit (EqO. Fitting with these six frequencies does not ade- 
quately describe the full variation of the star, but is sufficient 
to determine the zero points of different sites (Table|4]) to- 
gether with the amplitudes and phases of these particular 
six frequencies (Table[2]). It is these amplitudes and phases 
which are used in mode identification. Differential F-band 
photometry from all sites adjusted for zero-point differences 
are shown in Fig. [5] Also shown is the fitted Fourier curve 
using the six frequencies. From the amplitudes and phases of 
the fitted Fourier frequency components, we find amplitude 
ratios and phase differences listed in Table. [3] and shown in 
Fig.El 



Table 2. AmpUtudes, A (mmag), and phases 4> (radians) for 
the three frequencies of highest amplitude: /i = 5.66106, /2 = 
7.14895, /s = 7.77557 d"! and harmonics of /j. The epoch of 
phase zero is JD 2455700.000. 
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J3 


10.8 it 0.7 




4/1 


7.9 ±0.6 


2.7543 ± 0.0806 






R 


/i 


141.7 ±0.7 


2.2406 ± 0.0050 


2/i 


44.9 ± 0.7 


2.4660 ±0.0155 


h 


17.3 ±0.7 


-0.3476 ± 0.0407 


3/i 


13.6 ± 0.7 


2.8603 ±0.0516 


/3 


9.5 ±0.7 


-0.5988 ± 0.0746 


4/i 


6.2 ±0.7 


2.6144 ±0.1124 
/ 


fi 


113.0 ±0.6 


2.1942 ±0.0055 


2/1 


36.2 ± 0.6 


2.4433 ± 0.0169 


h 


13.7 ±0.6 


-0.4963 ± 0.0449 


3/1 


10.6 ±0.6 


2.8294 ± 0.0579 


/3 


8.3 ± 0.6 


-0.6936 ± 0.0746 


4/i 


5.3 ± 0.6 


2.6109 ±0.1147 



The harmonics of the principal frequency, fi , have large 
amplitudes and can be used as additional data for mode 
identification. In fact, the amplitude ratios and phase differ- 
ences for the harmonics 2/i , 3/i , are exactly the same as for 
the fundamental, /i, as can be seen in Fig.|3] 



5 MODE IDENTIFICATION 

There are several techniques that can be used to discrimi- 
nate between different values of the spherical harmonic de- 
gree, I. In the simplest technique, the amplitudes ratios rel- 
ative to a given waveband are compared with the calculated 
amplitudes ratios. The method relies on the fact that the 
variation of light amplitude in different wavelength bands 
depends on I and is independent of the azimuthal number 
m. It should be noted that our current models for photo- 
metric mode identification are only valid in the limit of no 
rotation. Rotation introduces departures from pure spheri- 
cal harmonics so that a single value of I is not sufficient to 
describe the eigenfunction perfectly. However, so long as the 
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Figure 4. Amplitude ratios (left panel) and phase difTerenees (radians, right panel) for /i and its harmonics, 2/i,3/i. 



rotation is not very rapid, a single value of I is still a reason- 
able approximation. We expect that for V2367 Cyg this will 
be the case, but the relative amplitudes will be modified to 
some degree due to this efTect. 

The most general expression for the magnitude varia- 
tion, mx, at wavelength A for a star pulsating with spherical 
harmonic degree, /, and angular frequency of pulsation, oi, 
with angle of inclination, i, is: 

Am;, = -i|^'^^r(cosi)6,;, {(1 - 0(/ + 2) cos(a;t) 

+ /T(aTA + Pt\) COs{ojt + 4jt) - fg{otgX + jSgx)} 

where e is the amplitude of the radial displacement, 

is the associated Legendre polynomial, /t = STcs/Tcs is 



the amplitude of local effective temperature variation, fg = 
5ge/ge is the amplitude of the local gravity variation for unit 
radial displacement at the photosphere and is the phase 
difTerence between effective temperature variation and radial 
displacement. The quantity bix = fihxPi{f^)dfi where the 
normalized limb darkening function is given by hx and /i — 
cos 8 is the polar angle in the spherical coordinate system 
where the z-axis is the line of sight to the observer. 

The first term within the bracket in the above expres- 
sion represents the change in magnitude caused by the geo- 
metrical variation of the star during pulsation and is inde- 
pendent of wavelength. The other two terms are wavelength 
dependent and lead to a variation of amplitude and phase 
with waveband. The difTerenees in amplitude and phase at 
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Figure 5. Light curves corrected for zero points showing fitted Fourier curve. 



different wavebands are small and accurate multicolour pho- 
tometry is required for good discrimination in I. Mode dis- 
crimination is improved by using wavebands covering as 
wide a wavelength range as possible. 

The quanti t ies oi t\ , Ptx, cigx, Pgx are defined in 
iBalona fc (| 19991 ) and can be determined from static 

model atmospheres. The assumption that at any instant the 
atmosphere of a pulsating star has the same structure as 
a static star of the same effective temperature and grav- 
ity, is probably a reasonable one if the frequency of pul- 
sation is low. At any rate this assumption is required in 
order to make progress with the current state of the art. 
The quantity, /t, the amplitude of local effective tempera- 
ture variation relative to radial displacement, is an impor- 
tant parameter in mode identification. It turns out that the 
calculation of /t is s ensitive to the treatm ent of convec- 
tion in & Set models (|Balona fc Everj[l999l ). One way of 
overcoming this difficulty is to leave /t as a free parameter 
to be determin ed by the observations using a m inimiza- 
tion technique (|Daszvnska-Daszkiewicz et al.l 12003 ). Appli- 
cation of this method to S Set stars indicates that the simple 
mixing-length theory for convection is inadequate, since the 



observed values of /t are not in agreement with the calcu- 
lated values. 

In mode identification we therefore have the option of 
leaving /t as a free parameter to be determined by the ob- 
servations themselves or to use the somewhat inadequate 
models to constrain the solution. We decided to use the lat- 
ter option because leaving /t free increases the number of 
parameters by two. We took the view that the effect of rota- 
tion, in any case, would limit the accuracy of the results so 
there is nothing to be gained by leaving /t as a free param- 
eter. The uncertainty in /t means, however, that not much 
importance can be placed in the resulting phase differences 
from the models. The amplitude of /t is more reliable and 
we can place more confidence in the calculated amplitude 
ratios. 

Balona et al. (2012) estimated the effective temperature 
Toft = 7300 ± 150 K and luminosity logL/Le = 1.7 ± 0.1, 
which places the star roughly in the middle of the instability 
strip and at the end of core hydrogen burning. There is a 
thin convection zone of H and He I near the surface which 
might modify the phase of the eigenfunction, but our current 
understanding of convection is insufficient to calculate this 
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Table 3. Amplitude ratios, A/ As, and phase differences (j> — (pB 
(radians) relative to the B band for six frequencies of highest 
ampUtude: /i = 5.66106, /2 = 7.14895, /s = 7.77557 d"!. 
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effect with any certainty. We may thus expect good agree- 
ment with the models for the relative amplitudes in various 
wavebands but less good agreement for the relative phases. 

The calculation of amplitude ratios and phase difTer- 
ences was performed using the FAMIAS software package 
l|Zimall2008l i. The photometric module of the software pack- 
age compares the observed parameters with a precomputed 
grid of non-adiabatic models for a range of stellar param- 
eters. The program produces amplitude ratios and phase 
differences in a number of different passbands which include 
the Geneva, Johnson/Cousins and Stromgren systems. In 
our case we used the results from the Johnson/Cousins sys- 
tem normalized to the B filter. 

Since the amplitude ratios depend to some extent on 
the stellar parameters, it is important to calculate these 
values over a sufficiently wide range of models centered on 
the best estimate of the parameters which we take to be: 
mass M/Mq = 2.10, effective temperature Tcff = 7300 K 
and \ogg — 3.5. We calculated amplitude ratios and phase 
differences for models in the range 7200 < Teg < 7400 K 
and 3.3 < logg < 3.7. These models were used because 



Table 4. Zero points in the UBVRI filters for various sites and 
the number of observations, A'^. 
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they produced the closest match to the observed frequen- 
cies /i,/2 and /a. In FAMIA S, models for 5 Set stars com- 
puted by the ATDH code JVentura et all l2008l l and MAD 
iMontalban fc Dupretl 120071 '). The results given by FAMIAS 
for different values of I are compared with the observed val- 
ues in Fig.O 

6 DISCUSSION 

At moderate rotation rates and for close frequencies, a 
strong coupling exists between modes with I differing by 2 
and of the same azimuthal order, m. A consequence of mode 
coupling is that modes of higher degree should be considered 
in photometric mode identification. Modes with nominal de- 
gree / > 2 acquire substantial I ^ 2 components and there- 
fore are more likely to reach detectable amplitudes. In these 
circumstances, mode identification using amplitude ratios 
or phase differences becomes both a spect- and m-dependent 
iDaszviiska-Daszkiewicz et al.ll2002l ). 

V2367 Cyg is a surprisingly rapid rotator. It can there- 
fore be expected that problems involving mode coupling will 
be significant and that non-adiabatic quantities calculated 
from non-rotating models will not be very reliable. In partic- 
ular, we expect that the observed amplitude ratios for /2 will 
be intermediate between I — and Z = 2 if mode coupling 
is inde ed present in this mode as suggested bv lBalona et al.l 
(|2012l ). 

From FigO we see that /i could be / = 0, 1 or 3. Unfor- 
tunately in this range of stellar parameters there is no dis- 
crimination between these three spherical harmonic degrees 
in the amplitude ratio. The phase differences are intermedi- 
ate between / = and 1 = 1. What we can conclude is that 
our results for /i are certainly consistent with I = and 
exclude 1 = 2. 

For /2 the amplitude ratios are intermediate between 
/ = and / = 2, as might be expected if this is a cou- 
pled mode. The phase differences are in agreement, though 
clearly no definite conclusion can be made on this basis 
alone. We conclude that mode coupling, as proposed by 
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iBalona et all \2Q1% . is consistent with our results. Results 
for /s are less reliable owing to the smaller amplitude, but 
are very similar to /2, so this might again be a coupled mode. 

Our overall conclusion from multicolour photometry is 
that mode identifi cation is in agreemen t with the identifica- 
tions proposed bv iBalona et al.l (|2012l ). but that the com- 
plexity introduced by the rapid rotation of V2367 Cyg pro- 
duces severe limitations on any interpretation. We certainly 
cannot find any evidence to contradict the identifications, 
which in itself is a useful result. 

Although HADS stars are of particular interest as they 
offer insights into nonlinear phenomena, the class itself is 
very poorly defined. There is no particular reason why the 
pulsation amplitude should be used as a criterion. Perhaps 
a better definition of the class is the presence of harmonics 
of the dominant frequency or of combination frequencies. It 
would also be very important to refine the physical parame- 
ters of HADS. If they are indeed intermediate between 5 Set 
stars and Cepheids, they should line in a region of the HR 
diagram between these two classes of stars. More precise ef- 
fective temperatures and luminosities are required to test 
this hypothesis. 
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